Black holes can have curly hair 
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We study equilibrium conditions between a static, spherically symmetric black hole and classical 
matter in terms of the radial pressure to density ratio p r /p = where u is the radial coordinate. 

It is shown that such an equilibrium is possible in two cases: (i) the well-known case w — > —1 as 
u — > Uh (the horizon), i.e., "vacuum" matter, for which p(uh) can be nonzero; (ii) w — » — l/(l+2fe) 
and p ~ (u — Uh) k as u — * Uh, where k > is a positive integer (w = —1/3 in the generic case 
k = 1). A non-interacting mixture of these two kinds of matter can also exist. The whole reasoning 
is local, hence the results do not depend on any global or asymptotic conditions. They mean, in 
particular, that a static black hole cannot live inside a star with nonnegative pressure and density. 
As an example, an exact solution for an isotropic fluid with w = —1/3 (that is, a fluid of disordered 
cosmic strings), with or without vacuum matter, is presented. 
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In real astro-physical conditions, black holes do not ex- 
ist in empty space but are rather surrounded by some 
kind of matter which is either in equilibrium with the 
black hole or is falling on it. In other words, real black 
holes are "dirty". Meanwhile, the famous no-hair theo- 
rems (see, e.g., P, HJ and references therein) are not di- 
rectly applicable to such situations of evident astrophys- 
ical interest. The main route in generalizing the possible 
black hole "hair" in such theorems consists in consider- 
ing different (dilaton, gauge etc.) fields, whereas a much 
simpler but physically and astrophysically more relevant 
environment, namely, macroscopic matter with certain 
pressure and density, drops out from consideration. 

The aim of the present paper is to partly fill this gap 
and to prove some statements of this kind in the frame- 
work of general relativity. Strange as it may seem, to the 
best of our knowledge, they were not found (or at least 
explictly formulated) before. 

The conditions we will rely on are the horizon regu- 
larity, the Einstein equations and the conservation law 
for matter. For simplicity, we restrict ourselves to static, 
spherically symmetric configurations. The manner of rea- 
soning is close to that of Ref. [|[ , where we have obtained 
some model-independent restrictions on the kinds of mat- 
ter able to support regular cosmological Killing horizons 
in Kantowski-Sachs geometries. 

We begin with writing the general static, spherically 
symmetric metric in the form 



where we have chosen the quasiglobal radial coordinate, 
corresponding to the "gauge" condition Opnoii = — 1- It 
has the following important properties [J, lf| : it always 
takes a finite value u = at a Killing horizon where 
A(u) = 0; 1 moreover, near a horizon, the increment 
u — Ufi is a multiple (with a nonzero constant factor) of 
the corresponding increments of manifestly well-behaved 
Kruskal-type null coordinates, used for analytic continu- 
ation of the metric across the horizon. Therefore, with 
this coordinate, the geometry can be considered jointly 
on both sides of a horizon in terms of a formally static 
metric (hence the name "quasiglobal"). On the other 
hand, both A(u) and r(u) should be analytic (or smooth 
at least up to derivatives of a certain order s > 2) func- 
tions of u at u — Ufa. A regular horizon corresponds to a 
regular zero of A(u), i.e., A(u) <~ (u— Uh) n , where n E N 
is the order of the horizon. In the case of a black hole, the 
outermost zero of A corresponds to the event horizon. 



Consider the Einstein equations 



—8wT" for the metric (JTJ, so that, due to the symmetry 
of the problem, the stress-energy tensor (SET) for an 
arbitrary kind of matter can be written as 

r ^ = diag(p, -p r , -p±,-p±), (2) 
where the density p, the radial pressure p r and the trans- 



ds 2 = A{u)dt 2 



du 2 

Mu) 



r 2 {u){de 2 + sin 2 9d(/> 2 ), (1) 
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1 In principle, u can take an infinite value at a candidate hori- 
zon where A — » 0, but then, as one can check, the canonical 
parameter of the geodesies also tends to infinity, so that the 
space-time is already geodesically complete and no continuation 
is required. Such cases, which can be termed "remote horizons", 
can be found, e.g., in some solutions of the Brans-Dicke theory 

We will not discuss them here. 

2 We use the units c = K = G = 1. 
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verse pressure p± are functions of u. 
One of the Einstein equations reads 

Gl-G\=2A r — = -^{p + p r ), (3) 

where the prime stands for d/du. Eq. ([3]) leads to a reg- 
ularity condition at u = Uh in terms of p and p r : 

Pr{u h )+p{u h ) =0, (4) 

or, more precisely, p r + p must have a zero of at least the 
same order as A(u) since, by regularity, \r"\ < oo. 

Regularity thus requires that, at a horizon, the null 
energy condition (NEC) 

> 0, = 0. (5) 

should be obeyed on the verge. Indeed, for the SET 
the NEC leads to 

Pr + P > 0, p± + p > 0. (6) 

The condition (TJ| is often discussed in connection with 
back-reaction of quantum fields (see, e.g., Q). Then, 
there is no great sense to speak of an equation of state; 
we simply have three different quantities in ob- 
tained from quantum mean values in a fixed background, 
as functions of the radial coordinate. In what follows, 
we discuss relations between them, mostly between p r 
and p, in quite a general form, and our reasoning is 
equally applicable to quantum mean values, classical 
fields and (which is astrophysically more relevant) usual 
matter with certain equations of state p r = p r (p) and 
P± =P±(p)- 

We are interested in general properties of matter sur- 
rounding a black hole horizon. A question of great im- 
portance for astrophysics is which kind of matter is con- 
sistent with the existence of a horizon. Although the 
condition ((4]) excludes p r and p being both positive, it 
can be satisfied if both quantities tend smoothly to zero 
as one approaches the horizon. This will be the subject 
of our study. Accordingly, mostly bearing in mind small 
p and p r , we will use the linear relation 

p r = wp, w = const. (7) 

We will assume p > 0. If, in addition, the NEC is 
satisfied (so that the weak energy condition is satisified 
as well) , we call the matter normal, otherwise it is said to 
be phantom. One can note that the NEC is often violated 
due to quantum effects [6(; phantom matter is also used 
in many studies as possible dark energy responsible for 
the accelerated expansion of the Universe. 

Consequences of the conservation law. The only 
nontrivial component of the conservation law V„T^ = 
can be written in the form (the prime denotes d/du) 

2r' A' 

Pr + — (Pr-p±) + ^-r(p + Pr)=0. (8) 



Suppose the validity of Eq. ([7]), at least near the horizon. 
As to the transverse pressure, we only assume that (at 
least, in the limit p — » 0) 

\p x \/p <oo. (9) 

It is a very weak restriction: indeed, for comparison, the 
dominant energy condition would require < 1. 

Then, near the horizon, the term with r' in Eq. ([8|) 
can be neglected as compared with the third one. In the 
leading approximation (i.e., retaining terms of the order 
p/Au), we obtain, as A — > 0, 

i4 -(ti>+i)/(a«) j w ^ . (10) 

The value w — (dust) is naturally excluded since non- 
interacting dust cannot be in equilibrium in a static grav- 
itational field. For different w, it follows from Eq. (fl"0|) : 

(i) w > (normal matter with p r > 0) or wj < — 1 
(phantom). The density diverges as A — > 0. Thus such 
matter cannot exist near a horizon. 

(ii) — 1 < w < (normal matter with p r < 0), in this 
case, both p and p r tend to zero at the horizon, and the 
condition (j4|) holds. 

(iii) w = —1: this special case requires more attention. 
First, Q now may hold with p r = —p =/= 0, which corre- 
sponds to a "vacuum fluid" considered in the next para- 
graph. Second, if we still assume p — > as A — > 0, one 
can check that such a solution to Eq. ([8]) near u — Uh can 
exist, but it cannot conform to the regularity require- 
ments connected with the Einstein equations to be dis- 
cussed below. Indeed, let us, going ahead, take p in the 
form (|15p and also assume a Taylor expansion of p r (p) at 
small p, 

Pr = -p + bp 2 + .... (11) 

Substituting all this into Eq. ([8]) and equating coefficients 
by equal powers of Au, we immediately obtain pk = 0, 
which means that p(u) cannot be represented by a Tay- 
lor series near u = Uh- It is a general result which 
does not depend on k or on the behavior of the func- 
tion r(p). One can find an asymptotic solution to Eq. (|8|) 
for small A under the assumption (fTTjl . having the form 
p w —2/ (bin A), so that really p — > as A — > 0; however, 
in accord with the above general result, this solution does 
not have the form (jT5j) and thus should be rejected. 

Inclusion of a vacuum fluid. Eqs. (|3]) and ([5]) for mat- 
ter under consideration do not change if we add a "vac- 
uum" matter with the SET 0] 

^(vac) = diag(p (vac ), Pr(vac) ; ~ RL(vac), _ 2>-L (vac) ) ) (I 2 ) 

and if there is no interaction between and T^( vac ) , that 
is, the conservation law holds for each of them separately: 
V„T£ = and V„T; (vac) = 0. 

The definitive property of vacuum matter is that the 
SET (|12p preserves its form under arbitrary radial boosts 
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Q. Examples of such matter are usual Maxwell ra- 
dial electric and magnetic fields for which pj_( va c) — 
P(vac)j their analogs in nonlinear electrodynamics with 
Lagrangians of the form L e = L e (F), F = F^F^ 0, 
Yang-Mills fields with a similar structure of the SET, 
and clouds of radially directed cosmic strings Q. In- 
dependently of a particular realization of such vacuum 
matter, a number of important properties follow from its 
algebraic structure, To( vac ) = T^( vac ) Q; it can be used 
both for constructing globally regular black hole models 
@, [|[ and for describing dark energy 0, [HI • 

By definition, vacuum matter does not contribute to 
the third term in ([8]). Therefore, the above conclusions 
are valid for matter with the SET ((2]) independently of 
whether or not there is a non-interacting admixture of a 
vacuum fluid with the SET (TT2")) , 

Consequences of the Einstein equations. There are 
two independent components of the Einstein equations 
for the metric (fT]), with the unknown functions A(u) and 
r(u). Assuming the total SET T£ (tot) = + T£ (vac) , 
we can choose such two components as Eq. @ and the 
equation containing only first-order derivatives, 



G\ = -^[-l + A'rr' + Ar' 2 } 



vac) ~Pr)- (13) 



Now, we require that the metric should be ana- 
lytic, or at least sufficiently smooth, in terms of the 
quasiglobal coordinate u (whose distinguished nature is 
discussed above, after Eq. (fT])) and thus admit continu- 
ation through the horizon. Therefore, we can write the 
Taylor expansions in Au = u — Uh'- 

A(u) = A n Au n [l + o(l% 

r{u) = r h +r' h Au+ \r'lAu 2 + o(Au 2 ), (14) 

with finite constants A n > 0, > 0, r' h and r'^. Recall 
that nSNis the order of the horizon. 

The l.h.s. of Eqs. ([3]) and (p~3j) are also smooth at the 
horizon. The same then applies to the r.h.s., in other 
words, both p (hence p r ) and /3( va c) are smooth and, in 
particular, since we are interested in configurations with 
p — > as u — > Uh, we can write in the same limit 



p = p k Au k [l + o(l)], 



k e N, 



(15) 



where pk — const > and k is the number of the first 
nonvanishing term of the Taylor series. Combining this 
with Eq. dTUJ) , we obtain: 



(16) 



2w 



2k 1 



: Vrl P\ 

-1 < w < 



This whole set 



a discrete set of values of w 
belongs to the range of interest, 

Let us now substitute these quantities to the Einstein 
equations. Eq. ([3]) in the main approximation (with the 
largest terms kept on each side) gives 



Evidently, finiteness of r'l leads to the requirement 



k > n, 



where k > n corresponds to r'l = 0. 



(18) 



Eq. (|13|> . in turn, gives in its main approximation 
[where the first terms on both sides are simply rewrit- 
ten while others are represented by their approximate 
expressions according to (TT4")) and (TTj 



- 1 + nA n r h r' h Au n 1 = -87rr^[p( vac) (w /l ) - wp k Au k ], 

(19) 

This relation leads to different results depending on the 
presence or absence of P( va c) • 



Indeed, if 



= at the horizon, the r.h.s. of (119 



2A n ^Au n 

Th 



-8tt(w + l)p k Au k 



(17) 



is zero at u — Uh, the only way to satisfy the equation 
is to require n = 1, and then A\rhr' h = 1. The horizon 
is simple, Schwarzschild-like. Evidently, the generic case 
in the expansion (fT5| is k = 1, and Eq. (JT6J) leads to 
w = —1/3 which, in case p r — p±, corresponds to a 
fluid of chaotically distributed cosmic strings (see [Tlj 
and references therein). 

In case /?( vac ) ^ 0, any n > 1 is admissible; if n > 1, 
Eq. CE!]) gives P( vac )(u h ) = l/(87rr£), while for n = 1 
we have P( vac )(u h ) = (1 - A 1 r h r' h )/(8irrl). Again, the 
generic case is certainly n = k = 1 and w = —1/3. One 
can also note that if matter is everywhere non-phantom, 
Eq. ([3]) leads to r" < 0, and if, in addition, there is a spa- 
tial asymptotic (not necessarily flat) r — -> oo as u — > oo, 
then r' > in the whole space, and r' h > in particular. 

We have actually proved the following theorem: 

Theorem 1. A static, spherically symmetric black hole 
can he in equilibrium with a static matter distribution 
with the SET f2]) only if near the event horizon (u — > Uh, 
where u is the quasiglobal radial coordinate) either (i) 
w — > — 1 (matter in this case has the form of a "vacuum 
fluid") or (ii) w — > —1/(1 + 2k), where w = p r j ' p and k is 
a positive integer. In case (i), the horizon can be of any 
order n, and p{uh) is nonzero. In case (ii), the horizon is 
simple, and p ~ (u — Uh) k . 

The generic case of such a non-vacuum hairy black 
hole is fc=l, implying w = — 1/3. In the isotropic case, 
Pr = P± ! it corresponds to a fluid of disordered cosmic 
strings Since such strings are, in general, arbitrarily 
curved and may be closed, one can express the meaning of 
the theorem by the words "non-vacuum black holes can 
have curly hair" . Recall, however, that in general our 
w characterizes the radial pressure, while the transverse 
one is only restricted by the condition ©. 

Other values of k (k = 2, 3 etc.) represent special 
cases obtainable by fine tuning of the parameter w. 

In the presence of vacuum matter with the SET (|12[) . 
the following theorem holds: 

Theorem 2. A static, spherically symmetric black hole 
can be in equilibrium with a non-interacting mixture of 
static non-vacuum matter with the SET {3P and vacuum 
matter with the SET ifT2j) only if, near the event horizon 
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(u — ► Uh), w = p r / p — > — n/ (n + 2fc), where u £ N is the 
order of the horizon, n < k € N, and p ~ (u — Uh) k ■ 

Thus a horizon of a static black hole can in gen- 
eral be surrounded by vacuum matter and matter with 
w = —1/3, which is true for any order of the horizon if 
n = k. (There also can be configurations with k > n and 
fine-tuned equations of state where w = —n/(n + 2k) > 
— 1/3.) An arbitrarily small amount of other kinds of 
matter, normal or phantom, added to such a configura- 
tion, should break its static character by simply falling 
onto the horizon or maybe even by destroying the black 
hole. In other words, black holes may be hairy, or "dirty" , 
but the possible kinds of hair are rather special in the 
near-horizon region: normal (with p r > 0) or phantom 
hair are completely excluded. In an equilibrium config- 
uration, all "dirt" is washed away from the near-horizon 
region, leaving there only vacuumlike or modestly exotic, 
probably "curly" hair. 

In particular, a static black hole cannot live inside a 
star of normal matter with nonnegative pressure unless 
there is an accretion region around the horizon or a layer 
of "string" and/or vacuum matter. 

We did not discuss the behavior of p± and p± ( vac ) (ex- 
cept for the restriction ©). In fact, these quantities 
are inessential for our reasoning but should be necessar- 
ily specified for finding complete solutions in particular 
models. Our inferences are quite general and hold for 
all kinds of hair: for instance, in all known examples of 
black holes with scalar fields (see, e.g., [l2| and refer- 
ences therein), the SETs near the horizon must satisfy 
the above conditions, which may be directly checked. 

Also, our approach is relevant to semiclassical black 
holes in equilibrium with their Hawking radiation (the 
Hartle-Hawking state), whose SET essentially differs 
from that of a perfect fluid. Since the density of quantum 
fields is, in general, nonzero at the horizon (see Sec. 11 
of the textbook [l[ for details), the regularity condition 
(|4]) tells us that such quantum radiation should behave 
near the horizon like a vacuum fluid. Our results show 
that a black hole can be in equilibrium with a mixture 
of Hawking radiation and some kinds of classical matter 
with — 1 < w < (including the important case of a Pas- 
cal perfect fluid with p r = p±). Possible effects of this 
circumstance for semiclassical black holes need a further 
study. Moreover, large enough black holes, for which the 
Hawking radiation may be neglected, can be in equilib- 
rium with classical matter alone, also including the case 
of a perfect fluid. 

Our reasoning was entirely local, restricted to the 
neighborhood of the horizon, and the results, which in- 
volve the single parameter w — p r /p, are in other re- 
spects mo del- independent. Meanwhile, a full analysis of 
specific systems would require the knowledge of the equa- 
tion of state and conditions on the metric in the whole 
space (e.g., the asymptotic flatness condition). Such an 
analysis depends on the model in an essential way and 
is beyond the scope of this paper. One can add that the 
equations of state well-behaved near the horizon are often 
incompatible with reasonable conditions at infinity (see 



the example below); it simply means that such matter 
does not extend to infinity and can only occupy a finite 
region around the horizon. 

It would be of interest to generalize our results to non- 
spherical and rotating distributions of matter. 

Example. In conclusion, let us present an exact so- 
lution for a system of utmost interest described by the 
above theorems. Consider a region of space with a 
non-interacting mixture of a vacuum fluid specified by 
87T/9( vac ) = A(it) [pj_(vac) is then found from the conserva- 
tion law for the SET l|12p] and an isotropic fluid of cosmic 
strings, such that p r = p± = —p/3- Then Eq. ([8]) leads 
to p = poA, Eq. © takes the form r" + a 2 r — 0, and, 
without loss of generality, we have 



r(u) = ro sin aw, 



(20) 



where pa > 0, ro > are arbitrary constants and a — 
(87rpo/3) 1/2 . The remaining unknown function A(u) can 
be found from Eq. (fT3"|) , which turns out to be linear, 



(21) 



A'rr' + A{r' 2 + a 2 r 2 ) = 1 - A(u)r 2 



hence easily integrable by quadratures for an arbitrary 
dependence A(u) (or A(r), as was used, e.g., in fft, [Io|). 
(Let us stress that our solution is different from that in 
Ref. @, obtained for a cloud of unidirectional strings.) 
In particular, for A = const we find 



A(u) 



1 — C cot au — Ar ( 1 — au cot au) 



, (22) 



with C = const. A horizon corresponds to A = 0, e.g., in 
case A = we obtain Uh — (1/a) arctanC; the horizon 
is simple: one can verify that A'(uh) ^ 0. As follows 
from ([2"0"]). this solution has no large r asymptotic, but 
it can be incorporated in an asymptotically flat model 
by matching it at some u > Uh to some intermediate 
layer (e.g., described by an analytic solution like the one 
for an incompressible fluid) admitting zero pressure at 
some surface, at which it can be further matched to the 
Schwarzschild solution. 

It seems instructive to trace the limiting transition 
from Eqs. (f^0|) . (|2"2"1) to the vacuum Schwarzschild-(anti) 
de Sitter metric. As the matter density vanishes, p — > 0, 
so that a — > 0, it is convenient to write ro = 0/ct, 
C = 2ma/j3, where and m are new constants. Then, 
after simple calculations, we obtain in the limit a — > 



r = flu, A(u) 



1 

J 2 



1 



2m 

u/3 



A 



-p 2 u 2 



Rescaling t i— > t — t/fl and using the coordinate r, we 
obtain the metric (1) with A(r) = 1 — 2m/r — Ar 2 /3, as 
required. 
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